Elliptical orbits in the Bloch sphere 
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As is well known, when an SU(2) operation acts on a two-level system, its Bloch vector rotates 
without change of magnitude. Considering a system composed of two two-level systems, it is proven 
that for a class of nonlocal interactions of the two subsystems including o~í ® Oj (with i,j 6 {x, y, z}) 
and the Heisenberg interaction, the geomètric description of the motion is particularly simple: each 
of the two Bloch vectors follows an elliptical orbit within the Bloch sphere. The utility of this 
result is demonstrated in two applications, the first of which bears on quantum control via quantum 
interfaces. By employing nommitary control operations, we extend the idea of controllability to a set 
of points which are not necessarily connected by unitary transformations. The second application 
shows how the orbit of the coherence vector can be used to assess the entangling power of Heisenberg 
exchange interaction. 

PACS numbers: 



I. INTRODUCTION 



The Bloch vector, or vector of coherence Q, provides a geomètric description of the density matrix of a spin-1/2 
particle which is commonly used in nuclear magnètic resonance. Mathematically, the Bloch vector may be viewed as 
. the adjoint representation of an su(2) object in an so(3) basis |2J. Extension of the notion of vector of coherence to 
çT} ' two-spin systems 

00, 

and more generally to quantum spin systems of higher dimensions |fj , has drawn attention in 
the contexts of quantum information theory and quantum computation. Specific motivations include the prospects 
of a uscful quantification of entanglement for composite systems 0, IS an d the quest for equations describing 
observables in quantum networks 0|. 

In the present work, the extension of the Bloch formalism to two spins is used to obtain a geomètric representation 
of the orbits of the vector of coherence for each spin system in the case that a nonlocal interaction of the form tji ® (Xj 
^Jj i is introduced. We propose that this geomètric picture will be useful in devising schemes for control of a quantum state 
via quantum interfaces i.e., through the mediation of an ancillary system. In this vein, we investigate the limits 
of control of a quantum state S, mixed or pure, given a nonlocal interaction and an ancilla Q. The simple geomètric 
picture developed below also applies to another special case of nonlocal interaction, namely the Heisenberg exchange 
Hamiltonian. As a second application of our formal results, we investigate the entanglement power of the Heisenberg 
interaction. 
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d • II. PRODUCT OF OPERATOR BASIS FOR A DENSITY MATRIX 



A. One qubit 



The density matrix p of a two-state system is a positive semi-definite Hermitian 2x2 matrix having unit trace. 
It can always be given expression in terms of the three trace-free Pauli matrices c^, i = 1, 2, 3, which are generators 
of su(2), and I/y/2 (I being the unit matrix): 

p=~I + va, (1) 

Here v is the vector of coherence, whose magnitude is bounded by <|| v ||< 1/2 because 1/2 < Tr(p 2 ) < 1. The 
two limiting vàlues of the norm correspond to maximally mixed and pure states, respectively. The magnitude of the 
Bloch vector differs by a factor of 1/2 from that of the vector of coherence, as a matter of convention. 

Unitary operations rotate the Bloch vector without changing its magnitude: SU (2) operations on the qubit 
correspond to 5*0(3) operations on the Bloch vector. The dynamical evolution of the Bloch vector under non-local 
operations is considered in the next section. 
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B. Two qubits and the correlation tensor 



In analogy to the representation Q), we adopt the generators of G = SU(4), i.e., the elements of the àlgebra 
g = su(4) (together with the unit matrix), as an orthonormal basis for the 4x4 density matrix of the two-qubit 
system. We employ this basis as it appears in Ref. 6, noting that it differs from the basis used in 0,0 only in the 
coefficients. 

The dynamical evolution of the system becomes more transparent if we choose basis elements of the àlgebra 
g = sit(4) in accordance with its Cartan Decomposition g = p(&e |fltll(|. The algebrasp and e satisfy the commutations 
relations 

[e,e]Ce, [p,e]Cp, [p,p]Ce. (2) 
The basis elements, Wj, j = 1, ... ,15 of the orthogonal àlgebra pair (e,p) are 

e = span-{cr x 1, a y ® 1, a z (g> 1, 1 ® a x , 1 ® a y , 1 ® a z } , (3) 



1 

p = span-{a x ® a x , <r x ® cTj,, a x ® o^, cr,, ® cr^, er^ <g) cr^, cr^ <g> a z , a z ® a x ,a z ® a yi a z ® a z } . (4) 
The basis defined by Eqs. © and |0J is used to expand the density matrix as 

15 15 

p = ^2Tr(pX j )X j = ^2p j X j , (5) 

where Xo = I/V4, Po = 1/V4, and Xj = —iWj (j = 1, . . . , 15). In this representation, the density matrix is specificd 
by three objects, namely the vectors of coherence i"i and r2 for the two subsystems along with the spin-spin correlation 
tensor TL 

í Pl \ í Pi \ í P7 PS P9 

ri = \ P2 , r 2 = p s , = I pia p n Pxi ) • (6) 

\P3 / \P6 ) \ Pl3 P14 Pl5 

We note that the object Tj has othcr names: Stokcs tensor ^ï|> entanglement tensor 0, and tensor of coherence 
(when combined with the coherence vectors in one object). Details of the properties of Tj can be found in Ref. |Ü, 
where many prior studies are cited. This tensor contains information on the correlations between the two subsystems, 
of both classical and quantum nature. Necessary and sufficient conditions for separability of a pure state can be stated 
in terms of its properties, whereas in the case of a mixed state, only necessary conditions for separability can be given 



III. EVOLUTION UNDER LOCAL AND NON LOCAL OPERATIONS 



As we have seen, the Lie àlgebra g — sm(4) possesses a Cartan decomposition g = e © p, which informs us that 
there exists within the Lie group G = S*[/(4) a subgroup of local operations Gl = SU{2) ® SU{2) generated by e. 
All the other operations are nonlocal and members of the coset space SU(A)/SU(2) SU(2), which does not form a 
subgroup of SU (4). It is known (see Proposition 1 of Ref. l·LOl·l that any U € S'ü" (4) can be written as 



U = k 1 Ak 2 



(7) 



with 



.4 



exp 



;( c i°Í cr * + c io\a 2 + c 3 al<7 2 z ) 



(8) 



where k\, k 2 £ SU(2) £g> SU(2) and ci, c 2 , C3 € R. In the following, we focus on the effect of nonlocal operations 
generated by a single operator among the possibilities for Ui ® where i,j 6 {x,y, z}. (This consideration includes 
the special case in which two of the parameters C\, c 2 , and C3 in the decomposition Q-© are zero.) Such nonlocal 
operations will be called one-dimensional. 
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(a) (b) 

FIG. 1: Local operations on the two spin subsystems produce a rotation of the corresponding vectors of coherence around some 
direction fi. The effect is the same for both pure states (a) and mixed states (b). 

A. Local operations 

Local operations are operations g E SU(2) ® SU(2) generated by the elements of e. From the commutation 
relations [e, e] C e and [p, e] C p it is clear that the elements of the vectors and tensor (T/) do not mix and do not 
affect one another. Under local operations, the vectors behave just like ordinary Cartesian vectors. In particular, a 
vector of coherence is rotated about some vector h as illustrated in Fig. 1, i.e., 

{r')\ = R)r j 1 , (r% = Ry 2 . (9) 

On the other hand, the correlation tensor transforms like a mixed Cartesian tensor, 

(T'Tj = R Ín( R ')" T n ■ (10) 

The magnitude of each object remains invariant under local operations. In addition, there exist fiftccn more invariants 
which can be constructed from the vectors and the tensor 121. 



B. One-dimensional nonlocal operations 

The nonlocal operations in the coset space SU(A)/ SU(2) ® SU(2) require, in their construction, exponentiation 
of at least one of the elements of p. Hence, under these operations the elements of the tensor and vectors of coherence 
are mixed, due to the commutation relations [p,e] C p and [p,p] C e. Wc shall establish that the one-dimensional 
nonlocal operations generated by the chosen interaction cr; ® Uj give rise to elliptical orbits for the vectors of coherence 
of the subsystems. The characteristics of these elliptic paths depend on the indices i and j, on the initial states of 
the subsystems, and on the degree of correlations between them. These orbits can be described by non-unitary 
transformations on each of the individual subsystems when one traces over the other's degrees of freedom. 

Accordingly, we take the interaction Hamiltonian between the two spins to be Hj — Ui® Oj/2, and, for reasons of 
simplicity, we suppose that the infernal Hamiltonians for the two spins may be ignored. Assuming that the duration 
of the interaction is <f>, and appealing to (i) the commutation relations as summarized in Ref. [10| and (ii) the identity 

exp [—i((j>/2)ai ® <7j] = cos((p/2)I — i sm(<p/2)<Ji <g) Oj , (11) 

we can make the following observations: 

(i) The components r\ and r 3 2 of the vectors of coherence remain unaffected; hence the vectors are confined to 
planes perpendicular to the i-axis and j-axis respectively. 

(ii) Of the nine elements of the correlation tensor TJ , only four experience changes. The five that are unchanged 
under the action of a <g> (jj are T„ ? and Tf with k ^ i and l ^ j. 
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(iii) The vectors r™ + Tj™ and r 2 ™ + Tf n are rotated, without change of magnitude, through an angle 4> about the i 
and j axes, respectively. (Here m ranges freely over {x, y, z}). 

(iv) More explicitly, the components of the vectors transform according to 



(r')\ = r{, 4^(r>%-, 2 . 



r 



í -> (r')i = rf cos - Tj sin , r f -> (r')™ = rg 1 cos - T^ sin , (12) 



(r')i = Tj° sin ^ + r[ cos , r£ -> (r') 2 l = T^ sin + r£ cos < 



and the components of the tensor of coherence, according to 

Tj -» (T')í = Tj , Tj — > (T')j = Tj , 

T/'^TOj-T/cos^-rjsin^, T^ - (T') l m = T^ cos - r% sin , (13) 
Tj -» (T')j = rf sin + Tj cos (j> , ï£ -> (T')j, = sin + T* cos tf> , 

with no change in the tensor's other elements. The ordered sets of indices (i,l,k) and (j,m,n) belong to 
{(x,y,z),(y,z,x),{z,x,y)}. 

Given this behavior, it is not difficult to show that ri(0) and r 2 (0)2 follow elliptical orbits. Since the 1,2 labeling is 
arbitrary, it suffices to demonstrate this property for the the vector ri (cj>). 

Proof. Referring to Fig. 2(a), the coordinates for a vector s tracing an ellipse in the x — y plane, with principal axes 
a and b rotated by an angle ip, are 

8x{<l>) = a> coscj) cosip + b sin0 sim/> , , , 

s y ((f>) = —a coscf> sxntp + b sin0 cosV> ■ 

The angle <fi is zero when the vector s is aligned with the principal axis a. 

The coordinates of the vector of coherence ri moving in the k — l plane are given by 

rf (0') = rf(0) cos 0'-Tj(O) sin 0', 

t\{tf) = T/(O)sin0' + r|(O)cos0'. 1 ' 

Of course, for the vector of coherence, 4>' = does not in general correspond to the principal axis a (see Fig. 2(a)). 
In fact, 4>' — <fi + Xi an d the coordinates of ri can be rewritten as follows in terms of the phase difference x : 

rf (0) = (rf (O)cosx- Tj (0) sinx) cos + (-Tj(O)cosx- rf (0) sinx) sin0, 

r[((j)) = (Tj= (0) sm x + r( (0) cx>s x) cos + (-r{ (0) sin x + T/(0) cos x) sin 0. ^ ' 

Comparison of the two sets of coordinates {s x ((/>), s y ((f>)} and {rf (<f>), r^ (0)} shows that a match can always be made, 
such that the parameters a, b, ip, and x can be determined by solving the system of equations 



a cos ip = rf (0) cos x — Tj (0) sin x , 

6sin^ = -Tj (0) cos x-rf(0) sinx, ^ 
a sin ?/> = — T fe (0) sinx — r^(0) cos x , 
6 cos ip = —r[ (0) sin x + Tj 5 (0) cos x ■ 



This completes the proof. It is important to note that the shape of the ellipse depends explicitly on the spin-spin 
correlation tensor. 

Solving Eqs. i|17|) for the angle x, we find 

, ,„ 2[rf(0)Tj(0)-r' 1 (0)Tf(0))] 

&n[ X ' -M(0)) 2 + (Tf (0))2 - (rf (0))* + (Tj(0))* ' [ 8j 

which specifies the initial orientation of the coherence vector ri with respect to the principal axis et. Suppose now the 
two-spin system is initially in a product state. For this case it is easy to prové these corollaries to our principal result: 

(1) The phase difference x is zero. This means that the initial positions of both coherence vectors lie on the a 
principal axis (as in Fig. 3(a)). It follows that the linear entropy of the state of each of the subsystems (defined 
by 1 — Trp 2 ) can only decrease, showing it is possible to increase the entanglement of the system with this 
interaction. (This will depend on initial conditions. See Section PVBI ) 
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(a) (b) 



FIG. 2: (a) The vector s describes an ellipse with principal axes a and b rotated by angle ip with respect to the y-axis. The angle 
4>, interpreted as the duration of a group operation, is measured relative to the a principal axis. (b) The vector of coherence r 
corresponding to one of the spins of the two-spin system moves on an ellipse on the k — l plane, with the angle (f)' measured 
relative to the original direction of r. 



(2) The length of the semi-minor axis of the ellipse followed by subsystem 1 is given by ò x = | r-^ (0)|[(y^(0)) 2 



(r^O)) 2 ] 1 / 2 (and likewise for subsystem 2 with 1 
pure state and the assumed single interaction <ii i 
vàlues of ir/2 and 3ir/2. 



-> 2 and {j, k,l} — > {i,n, m}). It follows that for an initially 
<7j , the maximum attainable entanglement is achieved at <p 



For the case of a initial state that is not pure but still separable, the phase difference \ does not vanish, in 
general (see figure 3(b)). Accordingly, the implied dynamical behavior of a classically correlated system distinguishes 
it from an uncorrelated system, but not from a system experiencing quantum entanglement. Moreover, the linear 
cntropy of each subsystem can either increase or decrease, showing it is possible to increase or decrease the amount 
of entanglement in the system. 

C. General nonlocal operations 

From Proposition 1 of Ref. Ilfl any nonlocal operation can be decomposed as a product of two local operations 
and an operation of the form 



A = exp 



(19) 



The operators {Yi} = i<TyUy/2, ia\a z /2} span a maximal Abelian subalgebra of P, and the relations 

1 „ 



-i\eijk\Yk 



(20) 



hold, where [•, •]+ denotes the anticommutator. Consequently, A of Eq. l|19fl can be written in product form, 



A = exp 



r(ci0- x er x ) 



exp 



exp 



:{c 3 a z a z ) 



(21) 
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(a) (b) 

FIG. 3: The initial position of the vector of coherence of subsystem 1 or 2 is shown, together with its time evolution under a 
one-dimensional nonlocal interaction (dashed line). If the initial state of the two-spin system is a product state, then the initial 
position is on the a principal axis of the elliptical path, as in (a). In general this agreement no longer occurs if the subsystems 
are initially correlated, either classically or quantum mechanically, as in (b). 



The property (|21|l telis us that any nonlocal operation can be decomposed into a sequence of operations effecting a 
succession of circular and elliptic paths in Bloch space. This result facilitates the calculation of the final states of the 
subsystems, but gives only limited insight into the geomètric characteristics of the coherence vectors' time orbits. For 
all Ci distinct, two general observations can be made: 

(1) The motion of the vectors of coherence is no longer restricted to a plane, since there is no linear combination of 
{a x gj 1, u y (8 1, o z ® 1} or {1 ® a x , 1 <g> a y , 1 (8 cr z } that is invariant under the action of A. 

(2) Characteristics of the trajectòries such as periodicity depcnd in detail on the parameters c\, c%, and C3. A 
trajectory is periòdic only if c^jcx and C3/C1 are both rational numbers. We note also that the set of parameters 
{ci, C2, C3} has been used to determine the equivalence classes of nonlocal interactions [ï(ï| as well as the invariants 
of the nonlocal interactions 0] . 



D. Special case of the Heisenberg Hamiltonian 

The Heisenberg exchange Hamiltonian, corresponding to c\ — c 2 = c 3 = —c/2, is not included in our general 
observations on nonlocal interactions (made for all Cj distinct), but like the one-dimensional Hamiltonians, it admits a 
simple geomètric picture. This interaction is the primary two-qubit interaction in several experimental proposals for 
quantum-dot qubits E^, E3> Il5| . It can also be used for universal quantum computing on encoded qubits of several 
types 0,0, 0,0,113, For these reasons, it warrants special attention. 

Introducing the time parameter <f>, the operator A of Eq. I|21|) now takes the form 

A{4>) = exp[—i(c(f>/2)(a x <g> a x + a y ® a y + a z ® a z )] 
= [cos 3 (c0/2) - isin 3 (c0/2)] I ® I 

-(i/2)e íc0/2 sin(c0)(o· x ® a x + a v ® a y + a z ® a z ). 

The time development of the density matrix under the operator A is given p(<fi) = A(<fi)p(Q)A* (<fi) and the corresponding 
coherence vectors change according to 

r\{4>) = |[rí(0) +4(0) + (rí(0) - 4(0)) cos(2c0) + (T^(0) - T*(0)) sin(2c0)] , (22) 
where k — 1, 2, 3 and cyclic permutations are implied. Similarly, for the coherence tensor we have 



= 2^(°) +^'(°) + ( T K°) - ï|'(O))cos(2c0) + (rf(0) - 4(0)) sin(2c0)] 



(23) 
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The elements of T2(4>) are found by symmetry 1 <-> 2. The quantities r\ + r 2 , Tj + T? , and T\ are unchanged by the 
operation, and the form of the one-parameter set that describes the time-evolving coherence vector is 



ri {4>) = R + S cos(2cç!>) + V sin(2cçi) , 



(24) 



where R = r x (0) + r 2 (0), S = n(0) - r 2 (0), and 



/T|(0)-T|(0) 
V= Tï(0)-T?(0) 



) 



(25) 



Clearly the vector traces out an ellipse lying in the plane spanned by S and V, defined by S x V, and passing through 
the point R. 



We shall now illustrate some of the results of Section III with two examples. The first provides a controllability 
result for nonlocal unitary interactions and the second demonstrates how the orbit of the cohcrcncc vector can be 
used to describe the entangling power of the Heisenberg exchange interaction. 



Let us now consider the implications of the findings of the preceding sections for the problem of quantum control. 
To this end, we adopt the nomenclature of Lloyd Q and identify spin 1 with the system S whose quantum state we 
wish to control, and spin 2 with the quantum controllcr or interface Q. It is assumed that (i) only one interaction 
Hamiltonian Hi is in play between S and Q and (i) system Q is complctcly controllable via control Hamiltonians 
{Hq} = {1 ® a x , 1 ® a y , 1 (g) a z } that span the su(2) àlgebra. The initial state of the bipartite system is taken to be 
a product state in the ensuing analysis. 

Suppose that the interaction Hamiltonian is nonlocal, but takes the one-dimensional form Hi = ai <Jj . Then 
the set {{Hq}, Hj} = {1 ® a x , 1 ® a y , 1 ® a Zl cr, ® a x , ai ® a y , ai ® a z } comprises a closed six-element subalgebra Gq 
of G. Given this set of operations, the vector of coherence rg of system S remains in the plane perpendicular to the 
i-axis. 

It has been established in Sec. III that when Hj — ai ® aj is the only element of the àlgebra su(4) affecting 
the two-spin system, the vectors of coherence ri and r 2 are constrained to move in elliptical orbits. Now, with the 
six-element subalgebra Gq available to the two-spin system S + Q, the reachablc set of the system S is enlarged to 
an elliptical disk (see Fig. 4). The principle axis of the disk coincides with the initial coherence vector rg(0) of the 
S system, while the length of its semiminor axis is given by b = [(^(O)) 2 + (rg(0)) 2 ] 1 / 2 |rg(0)|, where rç(0) is the 
initial coherence vector of system Q. 

Proof. First, if one implements the two-step sequence of a local operation € 1 <S> su(2) on system Q followed by the 
nonlocal operation Hj = a- t ® aj on S + Q, the orbit of r$ is necessarily an ellipse whose a principle axis lies along 
the initial coherence vector rg(0) and whose semimajor axis b is restricted by < b < [{rg(0)) 2 + (' , 5(0)) 2 ] 1 ^ 2 |rQ(0)|. 
Hence the set reachable by this two-step procedure is the elliptic disk in question. Second, using the Baker-Campbcll- 
Hausdorff formula one can show that all the elements of the six-element subalgebra G§ can be constructed by this 
two step sequence, so their reachable sets are the same. 

From this result we infer that the entropy of system S cannot be decreased by intervention of the quantum interface 
Q if the interaction Hamiltonian is limited to the form Hj = ai® aj. Noting that \yq\ < 1/2, it follows that a > b, 
where a and b are respectively the magnitudes of the semimajor and semiminor axes of the elliptical reachable set. 
Furthcrmore, it is secn that the systems S and Q become maximally entangled if the initial state of the system S is 
situated on the equatorial plane perpendicular to i-axis. 



Upon examining Eq. I|22l) . we see that the maximum entanglement, realized in a maximally entangled pure state, 
can be achieved if ri(0) = — r 2 (0), \ri\ = 1/2, and C(j) = ±7r/4. Otherwise, the state is not perfectly entangled since 
the linear entropy 1-Tr(p 2 ) is not minimized. This conclusion agrees with the result of Zhang et al. |10| that the only 
perfect entanglers that can be achieved with the Heisenberg Hamiltonian are the square-root of swap and its inverse. 



IV. APPLICATIONS 



A. Quantum control via quantum controllers and one-dimensional nonlocal interactions 



B. Entanglement power of Heisenberg interaction 



n Reachable set 



FIG. 4: The gray area is the set of reachable states for the system S if one has full control of the controller Q and the interaction 
Oi <g> Oj is available. This elliptical disk is characterized by a semimajor axis coincident with the initial vector of coherence for 
S and a semiminor axis with b = [(r|(0)) 2 + (r j, (0)) 2 ] 1/2 |r Q (0) | . 

However, suppose that the initial state of the two-spin system is represented by 

p(0) = ^(I + <r z )(g> (I + a z ), 

which is a pure-state density matrix for which r\ = 1/2 = rf and T* = 1/2, all other elements of the coherence 
vectors and coherence tensor being zero. Thcn 

íf (0) = 7^(0) cos 2 (c0) +r z 2 (0) sin 2 (c0) , 

while all other components of ri and r 2 vanish at time </>, and all other rf (0) = 0. In this case the ellipse collapses to a 
linc and the coherence vector simply oscillates between two vàlues along that line. The only element of the correlation 
tensor that changes is 

which vanishes for an initial tensor product of pure states for which the subsystems are polarized in the +z direction. 
Therefore one cannot create maximally entangled states with these initial conditions. 

V. CONCLUSIONS 

In this paper we have developed a geomètric representation for the orbits of the coherence vectors of a two-qubit 
system. In various circumstances we have shown that their evolution is described by elliptical orbits lying within the 
surface of the Bloch sphere. Importantly, every two-qubit unitary operation can be expressed as a combination of one 
of the evolutions we have considered, together with "pre" and "post" local single-qubit rotations. We anticipate that 
this geomètric picturc will be hclpful in devising schemes for control of a quantum state via quantum interfaces, and 
we have obtained a controllability result appropriate for such applications. Given the utility of the coherence- vector 
picture for modeling quantum systems and describing their entanglement, further studies along similar lines may be 
fruitful. Such work could include analysis of the orbits of higher-dimensional quantum states, as well as consideration 
of the effects of measurement operations on controllability. 
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